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Abstract. Despite the plethora of deep (sub-mJy) radio surveys there remains considerable doubt as to the exact nature of 
the galaxies contributing to the source counts. Current evidence suggests that starformation in moderately luminous normal 
galaxies is responsible for the bulk of the emission below lmJy. However given the sensitivities of these surveys we would 
expect a fraction of these sources to be distant radio galaxies. Using deep VLA and GMRT data we have found ~ 20 high-z 
candidate radio galaxies in two fields using the classical ultra-steep radio spectrum technique (De Breuck et al. 2000) and 
selecting galaxies with faint (i' > 25) optical counterparts. Several of these sources have X-ray detections in our deep 
XMM/Chandra observations and have fluxes high enough to put them in the quasar regime if they lie above redshift 3. 
Recently performed Spitzer GTO observations and upcoming near-infrared observations will help reveal the nature of these 
sources. 
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1. Introduction 

Deep extra-galactic radio surveys have been detecting numer- 
ous galaxies at 20cm with fluxes below lmJy for 20 years 
now (Condon et al. 1989; Windhorst et al. 1990; Hopkins et 
al. 1998; Ciliegi ef al. 1999; Gruppioni et al. 1999; Richards 
2000; Prandoni et al. 2001; Bondi et al. 2003; Hopkins et 
al. 2003; Seymour et al. 2004). However understanding the 
nature of these sources has not progressed far. Classically ra- 
dio surveys have detected quasars (and radio galaxies if the 
central black hole is hidden by dust in the optical) where the 
radio emission is often dominated by huge radio lobes now 
thought to be powered by jets. As surveys became more sen- 
sitive weak radio emission was found to be associated with 
starforming galaxies (eg Condon 1992) and that it was corre- 
lated with the far-infrared emission (Carilli & Yun, 2000). 

1.1. Source counts 

Source counts are the most basic diagnostic of any survey of 
the sky and give us our first chance to interpret the data we 
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have collected. The many surveys mentioned above all find an 
up-turn below lmJy in the Euclidean normalised differential 
source counts above that expected from the extrapolation of 
the counts at higher flux densities. This has been interpreted 
by many authors as being due to a new population of normal 
star forming galaxies and there is some evidence that starfor- 
mation is the dominant emission mechanism in the radio for 
several sub-mJy sources: 

1. Our two deep VLA 1.4GHz surveys (referred to as the 
13hr and Ihi XMM/RO SAT deep survey fields) have found 
~ 450 sources with flux densities greater than 30 /xJy in 
each 30' diameter region. We have found the number of 
unambiguous Active Galactic Nuclei, AGN, (ie large ra- 
dio lobes and/or broad optical emission lines) decreases 
regularly with decreasing flux density (Seymour et al. 
2005, in prep.). 

2. Modeling of the source counts finds that the starforming 
population of galaxies needed to explain the sub-mJy up- 
turn must evolve quickly with redshift (Seymour et al. 
2005). This evolution is compatible with (and constrained 
by) the need for the starformation rate density to rise 
quickly from z = — 1.0 (Hopkins et al. 2004). 



2 OUR GMRT DATA 
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Fig. 1. R— band distribution of the magnitude of optical coun- 
terparts to the our I3hi XMM/ROSAT 1.4GHz survey, includ- 
ing the 165/449 that have spectroscopically confirmed red- 
shifts. 

3. Alfonso et al. 2005 have found direct evidence from spec- 
tral line ratios for starformation rates high enough to ex- 
plain the observed radio luminosities of these faint galax- 
ies. They also find that luminosity function of the star- 
forming radio population evolves up to z ~ 0.5 at rate 
consistent with the models and analysis (e.g. Haarsma et 
al. 2000; Hopkins 2004) 

However the identifications of these radio sources are 
limited by their optical magnitudes. Figure ^ shows the dis- 
tribution of R— band magnitudes of optical counterparts to 
sources from the 13hr field. As can been seen there is a con- 
tinuous distribution sources from R ~ 18 — 25 before a tail 
off due to optical incompleteness. Around 20% of sources re- 
main un-identified, similar to other comparatively deep sur- 
veys. Also shown in Fig.^is the fraction sources with spec- 
troscopic identifications, less than half of the total sample. 
Other published surveys have similar spectroscopic detection 
rates. Therefore the current conclusions on the nature of the 
faint radio source population are drawn from only the opti- 
cally brighter sources. 

An additional issue is that some sources may have con- 
tribution to the radio from both AGN and starformation as 
can be seen in the local universe. In Fig. [2] we can see ra- 
dio contours overlayed onto an optical image of the Seyfert 
1 galaxy NGC405 1 . As well as the core being very luminous 
in the radio, there is also more extended radio emission from 
the starformation in the spiral arms. Although in this case the 
AGN outshines the radio emission from starformation by an 
order of magnitude, it does not mean that in other galaxies 
with less active black holes and weaker radio jets that the 
AGN and starformation will not be radiating more equally. 

2. Our GMRT data 

Radio spectra provide a particularly useful diagnostic. Star- 
burst emission is optically thin synchrotron emission, having 
a normal (i.e. steep, a ~ —0.7 where S oc v a ) spectrum, 
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Fig. 2. NGC 4051 at 20cm (contours) from a combination 
of VLA C & D Array overlaid on a Digital Sky Survey 
grayscale. It is clear that as well as strong radio emission from 
the AGN there is significant contribution from starformation 
in the spiral arms. 



and is also likely to be extended on galaxy-size scales. Pow- 
erful radio AGN, with extended steep spectrum radio lobes, 
are easily distinguished, but are rare in deep surveys. More 
commonly, emission from faint AGN is point-like and of rel- 
atively flat spectrum, perhaps from a compact jet. However 
some older AGN will have diffuse steep spectrum lobes that 
will be brighter at lower frequencies and hence be missed by 
high-resolution observations at 20cm. 

Therefore in august 2004 we obtained data on our lhr 
and 13hr fields at 610MHz and 235MHz with the Giant Me- 
trewave Radio Telescope (GMRT) in order to obtain 3-point 
radio spectral indices of our sources. From 4hours of ob- 
servations, including calibration, we have obtained a map 
at 610MHz with a rms noise level of 57/xJy/beam at the 
phase centre with a dynamic range just below 1000 making it 
the deepest 610MHz extra-galactic survey to the best of our 
knowledge (Moss et al., in prep.). 

Over 200 sources are detected in each field at 610MHz 
within a diameter of 1° with fluxes greater than 300 iiiy. 
Around 100 of these sources are also detected in the smaller 
1.4GHz VLA maps. Most steep spectrum sources, i.e. a ~ 
—0.7, are associated with optically bright spiral galaxies. We 
have several ultra-steep spectrum sources, a < — 1, with faint 
optical counterparts. The distribution of %'— band counterparts 
to these 610MHz sources is more bi-modal (see Fig.|3jl than 
the distribution of optical counterparts to the 1 .4GHz surveys. 
This distribution is consistent with the observations that many 
sources are bright spirals, but what is the nature of the opti- 
cally fainter sources? 



3.2 Expected Redshift - Luminosity Distribution of our sample 
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Fig. 3. i'— band distribution of the magnitude of optical coun- 
terparts to the lhr XMM/ROSAT 610MHz survey. 

3. HzRG candidates 

Both radio galaxies and quasars (i.e. radio-loud QSOs) host 
massive black holes and are very luminous at all wavelengths, 
but radio galaxies are preferred for the study of galaxy for- 
mation at high redshifts as quasars tend to out shine their 
host galaxies. In radio galaxies the active nucleus is hidden 
by optically thick material just beyond the accretion disk and 
so the host galaxy is detectable. There is now considerable 
evidence that high-z radio galaxies (HzRGs) are hosted by 
massive galaxies at every epoch (Archibald et al. 2001; Pen- 
tericci et al. 2001 ; Stern et al. 2003; Rocca-Volmerange et al. 
2004) 

High redshift radio galaxies are rare but, over the last 
decade or so, it has become clear that ultra-steep radio spec- 
tra (USS; a < —1.3; S oc v a ) provide the best selection 
diagnostic of high redshift radio galaxies (e.g. De Breuck et 
al. 2000). In the De Breuck et al, sample 35% of USS can- 
didates are confirmed to have z > 3. The De Breuck et al. 
sample has a flux limit of Siaghz > 10mJy. This sample 
preferentially samples high luminosities in the redshift range 
1 to 5. But what will we find at lower flux limits? 

There are two main possibilities. We may find lower lumi- 
nosity sources at high redshifts, and so be able to study prop- 
erly the evolution of the radio luminosity function at high red- 
shifts. We may also find HzRGs at unprecedented redshifts 
(ie z > 6) and so observe galaxies and massive black holes in 
their earliest stages of formation. Thus USS samples selected 
at low (sub-mJy) flux limits have great potential for the study 
of galaxy formation. 

3.1. Our Sample 

In each field we find around 10 sources with a < —1.0 and 
no confirmed optical counterparts, i.e. i' > 25. These make 
up our candidate sample which include several objects also 
detected in qmxXMM/ Chandra observations. If these sources 
are high-z radio galaxies then the X-ray emission is most 
likely from the accretion onto the central supermassive black 
hole. Whilst X-ray emission from the AGN is often absorbed 



10 



8 - 



>- 6 - 

0) 

E 

D 

z 4- 



2 - 



I 

-2.5 -2.0 -1.5 -1.0 -0.5 0.0 0.5 

Spectral Index 

Fig. 4. Distribution of radio spectral indices (al'iQMH Z , S u cx 
v a ) of sources common to the 1.4GHz VLA and 610MHz 
GMRT surveys of the lhr XMM/ROSAT field. 

in radio galaxies, mainly at harder energies > IkeV, the ab- 
sorption may be shifted out of the hard X-ray band at high 
redshifts, partly negating the decrease in flux density from 
the luminosity distance squared dimming. 

3.2. Expected Redshift - Luminosity Distribution of our 
sample 

The redshift-luminosity (1.4GHz) distribution of 165/449 ra- 
dio sources in our 13hr field, representing ~ 80% of the 
sources with R < 22, is presented in Fig.|5]Confirmed AGN 
(broad emission lines and/or radio lobes) are plotted as stars 
and the rest, presumed to be star-forming galaxies although 
there are probable exceptions, are represented by filled cir- 
cles. Sources circled have also been detected with our deep 
XMM/Chandra data (which is identical to that of the lhr 
field). Plotted as reference is a sample of high redshift ra- 
dio galaxies selected to cover the typical redshift-luminosity 
space of HzRGs and currently being studied with Spitzer in a 
GO Cycle 1 proposal (PI Daniel Stern). The detection limit of 
30/iJy is indicated by the dashed line. The expected range of 
any detected high-z radio galaxies is indicated by the shaded 
region. The boundary in redshift is determined by z = 3 cut- 
off due to USS selection and the faint optical nature of the 
sources, whilst the z = 7 cutoff represents the approximate 
era of reionisation and the earliest a supermassive black-hole 
can have gained enough mass to be radiating at the proposed 
luminosities. The boundary in luminosity is determined sim- 
ply as being 2 orders of magnitude from the detection limit 
at 1.4GHz ( given the pre-selection of Sqiqmhz 

> 0.3mJy), 

i.e. within the reasonable expectation of the value of the lu- 
minosities due to Malmquist bias. 

Figure [5] therefore shows that any high-z radio galaxies 
we detect will be either less luminous counterparts to the 
known HzRGs (although too luminous to be star-forming 
galaxies) or will be HzRGs at un-precedented redshifts likely 
accreting at high-rate as early super-massive black-holes are 
being formed. In the former case detecting less luminous 
counterparts to known HzRGs will increase the number of 
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Fig. 5. Distribution in luminosity /redshift space of the 165 ra- 
dio sources from our 13hr field with spectroscopic redshifts 
(stars = AGN, filled circles = starforming galaxies, circles = 
detection by Chandra and/or XMM). The sample of high-z 
radio galaxies is that currently being studied by Spitzer with 
GO Cycle 1 data (PI Daniel Stern) and is indicative of the 
typical luminosity distribution of HzRGs above z = 1. The 
shaded region shows the expected distribution of our candi- 
date sources given the select effects (radio flux and no optical 
ID) and the likely impact of the onset of reionisation at z = 7. 

low-luminosity HzRGs improving the study of observables 
influenced by the AGN power (e.g. the amount of scattered 
light). Additionally these observations would be vital in con- 
straining the less luminous end of the radio AGN luminos- 
ity function (Dunlop & Peacock 1990; Willot et al. 2001; 
Brookes et al. , these proceedings). In the latter case we may 
have detected spheroidal galaxies and super-massive black 
holes co-forming (e.g. Magorrian et al. 1998; Tremaine et al. 
2002). 

4. Future work & conclusions 

We have demonstrated the feasibility of using sub-mJy, multi- 
wavelength observations to select HzRG candidates using 
the VLA and GMRT These sources may either be low lu- 
minosity counterparts to the traditional family of HzRGs or 
HzRGs at unprecedented redshifts still in the early stages of 
formation. Either way these sources can provide vital clues 
to the early evolution of SMBH and giant elliptical galaxies. 
Clearly these sources warrant follow-up and through a com- 
bination of guaranteed UKIRT WFCAM and Palomar WIRC 
time we will obtain NIR imaging of these sources to K ~ 21 
in the coming semester. We have also applied for deeper 
GMRT time on each field to an rms sensitivity of 25 /^Jy at 
610 MHz in order to use radio spectra to help us distinguish 
between AGN and starforming galaxies. These deeper obser- 
vations will no doubt throw up more HzRG candidates. 

Additionally in June of this year the 13hr field was ob- 
served at Infrared wavelengths with the IRAC and MIPS in- 
struments aboard Spitzer as part of the instrument teams' 
GTO time (Pis Rieke and Fazio). These observations will 



be of an equivalent depth to the 4th deepest Spitzer survey. 
These data are clearly crucial in probing the restframe optical 
NIR emission of these potential HzRGs and determining ini- 
tial photometric redshifts. Spectroscopic redshifts could later 
be obtained through either NIR spectroscopy if bright enough 
or through the IRS instrument on Spitzer {e.g. Houck et al. 
2005). 
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